Freeze-out configuration in multifragmentation 
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The excitation energy and the nuclear density at the time of 
\ breakup are extracted for the a + 197 Au reaction at beam en- 

ergies of 1 and 3.6 GeV/nucleon. These quantities are calculated 
q ■ from the average relative velocity of intermediate mass fragments 

(IMF) at large correlation angles as a function of the multiplic- 
ity of IMFs using a statistical model coupled with many-body 
^ \ Coulomb trajectory calculations. The Coulomb component v c 

i— h \ and thermal component Vq are found to depend oppositely on 

1 the excitation energy, IMFs multiplicity, and freeze-out density. 

Ph ■ These dependencies allow the determination of both the volume 

>- . and the mean excitation energy at the time of breakup. It is 

found that the volume remained constant as the beam energy 
| was increased, with a breakup density of about po/7, but that 

the excitation energy increased 25% to about 5.5 MeV/nucleon. 
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In spite of intensive experimental and theoretical efforts, the mechanism 
of nuclear multifragmentation remains a subject of much interest and debate. 
The study becomes even more attractive with the availability of more high 
quality data from powerful 4ti detectors [|TJ. 

In this Letter we report on results of the determination of the freeze-out 
density and excitation energy by studying the average relative velocity 
< v rei > of IMFs at large correlation angles as a function of the multiplicity 
of IMFs. We use the Berlin version of a statistical model || coupled with the 
many-body Coulomb trajectory calculations to describe the multifragmenta- 
tion process. Among the many available multifragmentation models, those 
based on the simultaneous breakup mechanism have been rather successful 
in explaining a number of aspects of multifragmentation data [[!]. Moreover, 
statistical phase-space models like ours have a unique ability for calculating 
many-fragment correlation functions with reliable statistics. 

In the simultaneous breakup scenario, one can picture the multifragmen- 
tation as follows: A composite nuclear system formed initially in any reaction 
will expand under thermal and/or compressional pressure. The expansion 
may take the system into the hydrodynamical instability region where the 
sound velocity is imaginary[Q. In the instability region the system breaks up 
into fragments and nucleons. At the same time the further expansion slows 
down due to the creation of new surfaces of fragments and becomes chaotic 
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due to the strong interactions among fragments and nucleons through the ex- 
change of nucleons and momenta. At the end of this chaotic expansion (the 
instant of freeze-out) the fragments have rearranged themselves such that the 
entropy of the system has reached its maximum. At the moment of the freeze 
out, the momentum distribution, spatial distribution, mass distribution and 
internal excitations of fragments of the system are solely determined micro- 
canonically by the available phase space. After the freeze out, fragments and 
nucleons move along their Coulomb trajectories. This picture is in agree- 
ment with a recent BUU-model simulation of the reaction dynamics leading 
to multifragmentation|5|. |J. 

Within the above picture, the most important quantities determining the 
final state of the multifragmentation events are the excitation energy and 
the density or size of the system at the freeze out. Ideally, these quantities 
can be calculated from dynamical models for nuclear reactions. However, the 
breakup conditions (time scale, entropy and partitioning into fragments) is 
the most interesting topic of current research because none of the existing 
dynamical models can handle the breakup properly. The extracted freeze- 
out density and excitation energy will put strong constraints on the reaction 
dynamics leading to the formation of hot nuclei and on the nuclear equation 
of state for low-density nuclear matter. Moreover, to a significant extent, our 
study itself is a study of the breakup dynamics. 
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The velocities of the fragments in the final state are sensitive to their spa- 
tial distribution at the freeze out due to their long range Coulomb repulsion. 
The velocities are also sensitive to the excitation energy of the fragmenting 
system. The charge and mass distributions of fragments, as well as the ther- 
mal velocities of the fragments at the freeze out, are mainly determined by the 
excitation energy of the system. Therefore, the relative velocity distributions 
are useful tools for studying the breakup process as well as for extracting the 
freeze-out density and excitation energy. The above discussions show that 
the two quantities have to be extracted simultaneously. It is interesting to 
mention that the relative velocity distributions for various fragment classes 
have been used to distinguish between the simultaneous breakup scenario 
and the sequential binary decay scenario []7|, ||, ||, |10| with contradictory con- 
clusions. By selecting dedicated event classes a deeper insight can be gained. 
For three-body events from the system Ar + Au at 30 MeV/nucleon it was 
shown that the relative velocity distribution can be well reproduced assum- 
ing the sequential binary decay of a hot system at an excitation energy of 
about 3 MeV/nucleon [0. This result can be understood as a deep inelastic 



collision, followed by the fission of the heavy partner [12| . At higher incident 
energies, selecting events involving only three heavy fragments, the ALADIN 
collaboration has very recently also studied the relative velocity distribution 
from 600 MeV/nucleon Au induced reactions [p}|, their results suggest a fast 
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disintegration of the entire highly excited system. A range of the freeze-out 
radii and excitation energies was extracted simultaneously using a simple 
fragmentation model coupled with three-body trajectory calculations. 

As a first part of this Letter we study the model predictions for the 
average relative velocity between IMFs with large correlation angles as a 
function of the total IMF multiplicity, varying the freeze-out density and 
the excitation energy. In the second part we compare these results with 
experimental data in order to deduce the freeze-out condition. This method 
goes far beyond the usual way of comparing the relative velocity distribution 
for events involving only three fragments or the relative velocity distribution 
obtained without any constraint on the multiplicity of IMFs. It is evident 
that the multiplicity and charge of the unobserved fragments have decisive 
influence on the observed velocity distributions. We want to emphasize that 
this means we test more than just a three-body dynamics as in ref . [|T3|] . 
The change of < v re i > with the IMF multiplicity at a relatively constant 
excitation energy gives the clue to our method of extracting the source size 
and excitation energy. This is a clear advantage using the very asymmetric 
collision system a + Au. There is also a qualitative difference between highly 
energetic a induced fragmentations and heavy-ion induced ones. In the later 
the variation of the IMF multiplicity is mainly due to the broad variation of 
the excitation energy. Here it is much more due to the fluctuation of Mjmf 
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at a relatively narrow distribution of excitation energies. 

The final velocity of each fragment v} can be decomposed as 

v} = vo + v c , 

where the Vq term is the thermal velocity of the fragment at the freeze out. 
At a fixed excitation energy of the fragmenting system the v term increases 
with the multiplicity of IMFs due to the smooth decrease of the average 
mass of the IMFs. The v c term is the velocity of the fragment gained after 
moving in the Coulomb field of the rest of the system and thus depends on 
the radius of the system at the freeze out. The dependence of v c on the 
multiplicity of IMFs is somewhat more complicated. It is strongly influenced 
by the existence of a third, larger fragment. This will be discussed in detail 
later. Therefore, the dependence of < v re i > on the multiplicity of IMFs 
results from the interplay between the two velocity components of the IMFs 
just mentioned. They reflect the freeze-out density and the excitation energy. 
The secondary deexcitation of primary hot fragments will slightly alter the 
relative velocity distribution. Since the deexcitation occurs late, the influence 
of the change in the Coulomb field is small. Further, since the emission is 
assumed to be isotropic, only the width will be changed, but not the mean 
value. 

Necessary inputs to our calculations are the mass and the charge of the 
fragmenting source. Since this Letter relates to a + Au reactions at 1 and 
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3.6 GeV/nucleon, cascade model calculations |14] were performed for these 



reactions, yielding the mass and charge of the residual nucleus, which is the 
emitting equilibrized source. For the higher incident energy a mean residue 
mass of A=150 and charge of Z=62 was deduced, for the lower incident en- 
ergy the mean residue mass and charge was determined to be A=170 and 
Z=70, respectively. The mass and charge of the heavy residue are rather in- 
dependent of the impact parameter due to the large mass asymmetry. This 
legitimates the comparison of the experimental results to the statistical cal- 
culations using a fixed residue mass, charge and excitation energy. This is 
an important advantage of the highly asymmetric a + Au system for inves- 
tigations of this kind. 

To see how the dependence of < v re \ > on the multiplicity of IMFs is 
sensitive to the breakup conditions of hot nuclei towards multifragmentation, 
we study in Fig. 1 (a) the excitation energy dependence of < v re i > at large 
correlation angles for a given breakup density. 

In Fig. 1 (a)-(c), the calculations are done for excitation energies of 4, 
5, 8 and 10 MeV/nucleon respectively at a fixed freeze-out radius of 2.1A 1 / 3 
and a residue mass and charge of A=150 and Z=62. It is seen that for 
small excitation energies < v re \ > decreases with increasing multiplicity of 
the intermediate mass fragments Mjmf and that < v re i > is higher at E* =4 
MeV/nucleon than at 5 MeV/nucleon. However, an opposite behaviour is 
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seen at high excitation energies. Now < v re \ > increases with increasing 
M IMF and < v ret > is higher at 10 MeV/nucleon than at 8 MeV/nucleon. 
These features of < v rd > can be understood by considering the interplay 
between the two velocity components of the fragments and their excitation 
energy dependence. For this purpose we show in Fig. 1 (b) the average 
maximum charge of the fragments as a function of the multiplicity of IMFs at 
the corresponding excitation energies. The dashed line at Z=15 is the upper 
boundary of the charge defined for the IMFs. At lower excitation energies 
(e.g. 4 and 5 MeV/nucleon) the average maximum charge is larger than 15 
for a significant range (if not all) of the IMF multiplicities. The fragment 
with the highest charge in each event accelerates the other fragments and is 
therefore responsible for the higher values of < v re i >. The maximum charge 
of these events decreases strongly with increasing IMF multiplicity. This 
explains the trend of decreasing < v re i > with increasing M IMF observed at 
lower excitation energies. 

At high excitation energies the maximum charge is within the charge 
range of IMFs for all statistically significant events (seen from Fig. 1 (c)) 
and the thermal term v dominates over the Coulomb term v c . Since the 
average mass of IMFs will decrease as the multiplicity of IMFs increases, 
the vo term and consequently < v re i > of IMFs will increase with increasing 
Mimf- The fact that < v re i > is higher at higher excitation energies is 
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due to both the higher temperature reached and the lower average mass of 
the fragments. To evaluate the statistical significance of each multiplicity 
of IMFs, we have presented the distribution of the multiplicity of IMFs at 
the corresponding excitation energies in Fig. 1 (c). The typical rise and fall 
behaviour of the IMF multiplicities as a function of the excitation energy can 
be seen between E* = 5 and 10 MeV/nucleon. The average IMF multiplicity 
at E* = 10 MeV/nucleon is lower than that at 8 MeV/nucleon. This is due 
to the increasing dissociation into light fragments (Z < 3) and nucleons at 
higher excitation energies. 

To see the effects of the freeze-out density or radius on < v re i > we show in 
Fig. 1 (d) < v re i > at large correlation angles as a function of Mimf for freeze- 
out radii of 2.1 A 1/3 , 2 AA 1/3 and 2.QA 1/3 respectively at a fixed excitation 
energy of 5 MeV/nucleon. The magnitude of < v re \ > increases rapidly 
with decreasing freeze-out radius due to the stronger Coulomb repulsion at 
smaller freeze-out radii. It is interesting to note that at R c = 2.6A 1 / 3 the 
Coulomb repulsion is so weak that for Mimf > 6 the thermal term v starts 
dominating over the Coulomb term v c and < v re i > starts increasing with 
increasing M IMF . 

After establishing the various dependencies of < v re i > we can compare 
them with the experimental data. For this purpose the reaction a + Au 
is ideally suited, since all fragments originate from one source, the target 
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residue. Additionally the recoil velocities are small which facilates the deter- 
mination of the velocity and angular correlations with high precision. The 
reaction has been studied at 1 and 3.6 GeV/nucleon incident energy at the 
synchrophasotron of the JINR in Dubna using the new 4n setup "FASA" [|l5|] . 
In |16|] it has been demonstrated that sufficient excitation energy is reached 
to induce the multifragment breakup of the system into many IMFs. 

The multiplicity of the IMFs was measured by a fragment multiplicity 
detector (FMD) system covering a large part of the 4ir solid angle. The 
energy, velocity and mass of single fragments from the event were deter- 
mined with high precision using time-of-flight telescopes (TOF). The relative 
angle and velocity of coincident fragments were measured using a position- 
sensitive parallel-plate avalanche counter (PPAC). For the FMD and the 
PPAC an IMF is defined as a fragment with a charge number of 3 < Z < 15. 
For TOF it is defined correspondingly as a fragment with a mass num- 
ber of 6 < A < 30. The large correlation- angle data were taken by 
three TOF-PPAC combinations covering the correlation angular range of 
105° - 155°, 130° - 180° and 150° - 180° respectively. Efficiencies of both the 
TOFs and PPACs are taken into account in the data analysis and the model 
calculations. 

The comparison to the experimental data is shown in Fig. 2, which gives 
the average relative velocity of the IMFs at large relative angles, detected 
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in TOF and PPAC, as a function of the multiplicity of the IMFs measured 
in the FMD. In this figure the efficiency of the FMD was not corrected but 
considered in the respective calculations. Figure 2 shows, that a combination 
of a freeze-out density of po/7 and an excitation energy of 4.5 MeV/nucleon 
for the lower incident energy and the same freeze-out density but a higher 
excitation energy of 5.5 MeV/nucleon for the higher incident energy can re- 
produce the experimental data reasonably well. With this set of parameters, 
our calculations also well reproduce many other aspects of the experiment 
simultaneously, including the average IMF multiplicity, the width and the 
shape of the relative velocity distribution as well as the mass distribution of 
fragments. To estimate fluctuations of < v re \ > originating from fluctuations 
of the residue charge and excitation energy we have performed calculations 
for the case of E beam /A = 3.6 GeV with Z=62, 70 and 79 respectively. By 
adjusting the excitation energy between 5.5 and 6.6 MeV/nucleon in order 
to have a similar average IMF multiplicitiy of about 5.9, the < v re i > varies 
by about 4% which corresponds to a variation of R c by about 5%, yielding 
an uncertainty of the breakup density of 15%. 

In summary, using a statistical model coupled with many-body Coulomb 
trajectory calculations we have studied the average relative velocity of inter- 
mediate mass fragments as a function of the multiplicity of IMFs at differ- 
ent but constant excitation energies and freeze-out radii of the fragmenting 
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source. These two quantities can be determined individually because of the 
complex and different forms of the dependence of < v Td > on the IMF 
multiplicity at different excitation energies and freeze-out densities. For the 
a+ 197 Au system at beam energies of 1 and 3.6 GeV/nucleon, the freeze-out 
density was found to be p c rs Po/7. While the mean excitation energy in- 
creases from 4.5 MeV/nucleon to 5.5 MeV/nucleon. The freeze-out density 
found is very similar to the ones deduced from nucleus-nucleus collisions. 
Conceivable compression effects in the later seem therefore to play a minor 
role if at all. 
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Figure Captions 

Fig. 1 (a) Average relative velocity of IMFs at large correlation angles 
as a function of the multiplicity of IMFs, (b) the average maximum 
charge of fragments as a function of the multiplicity of IMFs, (c) the 
multiplicity distribution of IMFs, and (d) average relative velocity as 
function of the multiplicity of IMFs for different freeze-out radii. 

Fig. 2 Dependence of the average relative velocity of IMFs at large correla- 
tion angles on the multiplicity of IMFs measured in the FMD (not cor- 
rected for efficiency) for the system a + Au at 1.0 GeV/nucleon (filled 
squares: experiment, solid line: calculation with E* / A = 4.5MeV and 
R c = 2.3* A 1 / 3 ) and 3.6 GeV/nucleon (filled circles: experiment, dashed 
line: calculation with E*/A = 5.5MeV and R c = 2.3 * A 1/3 ) respec- 
tively. 
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